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Abstract 

In the precedent paper of the authors (hep-ph/0510410), the bb states were treated in 
the framework of the spectral integral equation, together with simultaneous calculations 
of radiative decays of the considered bottomonia. In the present paper, such a study is 
carried out for the charmonium (cc) states. We reconstruct the interaction in the cc- 
sector on the basis of data for the charmonium levels with J-^*^ = ^, 1 , 0^+, l^"*", 
2++, 1+^ and radiative transitions ip{2S) ^Xcoi'^P), 7Xci(l-P), jXc2{^P), "yVci^S) and 
Xco(l-P)) Xcii^P), Xc2(l-P) — > 'jJ/ijj. The cc levels and their wave functions are calculated 
for the radial excitations with n < 6. Also, we determine the cc component of the photon 
wave function using the e"'"e~ annihilation data: e"'"e~ J/'i/'(3097), VlSOSG), -0(3770), 
■0(4040), ■0(4160), ^0(4415) and perform the calculations of the partial widths of the two- 
photon decays for the n = 1 states: r]co{lS), Xcoi^P), Xc2(l-P) — ^ 77, and n = 2 states: 
rico{2S) — > 77, Xco{2P), Xc2(2-P) — > 77. We discuss the status of the recently observed cc 
states A(3872) and y(3941): according to our results, the A(3872) can be either Xcii^P) 
or ric2{lD), while y(3941) is Xc2(2P). 



1 Introduction 

In this paper, we continue the study initiated in [1] for the bb states in the framework of the 
spectral integral equation. Here the results for the cc-states are presented. 

In [2], the program has been formulated for the reconstruction of the quark-ant iquark 
interaction based on our knowledge of meson levels and their radiative decays. Within this 
program, as a first step we have considered the bottomonia [1]. Now we present analogous 
results for the charmonia. In the subsequent publication, we suggest to give corresponding 
results for the light quark-ant iquark systems. 
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Our study is carried out in terms of the spectral integral technique. The application of 
this technique to the composite quark-antiquark systems and its relation to the dispersion 
N / D-method have been discussed in [1,2] — one may find there necessary details. 

Still, let us point once again to particular properties of our approach. The quark-antiquark 
interaction given, the spectral integral equations provides us unambiguously with both levels 
and wave functions of composite systems. But if the interaction is unknown, to reconstruct it 
one needs to know the levels as well as their wave functions. Our knowledge of the interaction 
of constituent quarks (in particular, long-range interaction) is rather fragmentary, so actually 
the description of the composite quark-antiquark systems means the reconstruction of quark 
interaction. To this aim, one needs the information on the wave function of a level, and the 
radiative decays are precisely the source of information for the wave functions. Because of that, 
in our investigation of the quark-antiquark states we rely equally upon our knowledge of levels 
and radiative decays. 

The method of calculation of the radiative transition amplitudes in terms of the double 
dispersion integrals was developed in a number of papers [3, 4, 5], where the important point 
was the representation of the transition amplitude in the form convenient for simultaneous 
fitting to the spectral integral equation — it was done in [6, 7]. 

An important information on the cc-meson wave function is hidden in the two-photon meson 
decays: cc — meson 77. For the calculation of such a type processes, the quark wave function 
of the photon is needed; correspondingly, the method of reconstruction of the 7 cc vertices 
was developed in [8]. The results of calculations of the cc-mesons are presented in Section 2. 
The known levels with J^^ = 0"+, 1~, 0++, 2++, 1+" have been included in the fitting 
procedure, together with the widths of radiative transitions as follows: ip{2S) — > ■jXcoi^P), 
7Xci(l-P), 7Xc2(l-P), -fr]c(lS) and Xco(l-P), Xci(l-P), Xc2(l-P) ^J/ijj. The masses of the states 
have been calculated for the radial quantum numbers n < 6 and radiative transition amplitudes 
for the states with n < 2. Let us emphasize that we face rather significant relativistic effects 
in the case of excited states. 

The recently observed states X(3872) and y(3941) [9, 10, 11, 12] are lively discussed at 
present time [13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. Our calculations argue that X(3872) is 
either the excited l+^-state Xci(2P), or basic 2~+-state ric2{lD). The charmonium y(3941) 
state can be the radial excited 2++ state, Xc2(2P). 

We determine the cc component of the photon wave function using the transitions e+e~ 
J/^/;(3097), V(3686), V'(3770), ^(4040), -0(4160), -0(4415) that allow us to calculate partial 
widths of the two-photon decays. We present the results for the IS*, IP states: ?7co(2979), 
Xco(3415), Xc2(3556) ^ 77, and 2S, 2P states: ?7co(3594), Xco(3849), Xc2(3950) ^ 77- The 
predictions are also given for the two-photon decays r^co, Xco, Xc2 77 of the states below 4 
GeV. 

In Conclusion, we briefly summarize the results. 
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2 Charmonium states 



We calculate the cc levels and their wave functions using two types of the t-channel exchanges: 
scalar and vector states, (/ /) and (7^ (8) 7/^). The calculations of the cc-systems have been 
carried out similarly to the consideration of bottomonia [1], with both retardation interaction 
(Solution i?(cc)) and three variants of instantaneous interactions (Solutions /(cc), II{cc) and 
U{cc)). In the case of instantaneous interactions, the t-channel exchanges may be represented 
by using the potentials. In the fitting to the cc spectra, we have applied scalar and vector 
potentials of the type: 

a + br + ce-''^'- + . (1) 

r 

The presentation of these interactions in the momentum space, also with retardation effects, is 
given in [1]. 

The interaction parameters obtained in the fit are as follows (in GeV units): 

Interaction Solution a b c /ic d /id 

0.172 1.141 0.201 0.000 0.401 

0.172 1.597 0.201 0.000 0.401 

0.150 -0.044 0.351 -0.245 0.201 

0.157 0.952 0.201 0.000 0.401 (2) 

0.000 1.225 0.201 0.506 0.001 

0.000 1.587 0.201 0.445 0.001 

-0.150 -1.600 0.201 0.544 0.001 

0.000 0.977 0.201 0.492 0.001 

As concerns the solution U{cc), we have included into calculations the scalar and vector con- 
finement forces, with bs = —by = 0.150 GeV^. The vector confinement potential is needed for 
the description of the light quark states (qq) with large masses. In this way, solution U{cc) 
gives us the description of data with a universal confinement potential for all fiavours. 

The as coupling, being determined by the one-gluon exchange forces, is of the same order 
in all solutions: = 3/4 • rfy ~ 0.38. 

The mass of the constituent c-quark is taken to be rric = 1.25 GeV. This mass value is 
consistent with the magnitude provided by the heavy-quark effective theory [23, 24]: 1.0 < 
rric < 1.4 GeV; a shghtly larger interval for rric is given by lattice calculations, 0.93 < rric < 1.59 
GeV, see [24] and references therein. The compilation [25] gives us 1.15 < rric < 1-35 GeV. 



(101) 



/(cc) 
//(cc) 

U{cc) 
R{cc) 



-1.980 
-1.980 
-0.300 
-1.551 



(7/. ® 7m) 



/(cc) -1.682 

//(cc) -1.682 

U{cc) 1.000 

Rice) -1.312 
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2.1 Masses of cc states 

The fitting procedure results in the following masses (in GeV units) for 1 states: 



State 


Data 


I(cc) 




nice) 




U(cc) 




Rice) 




Rj 


J/ih 


3.097 


3.100 


(S) 


3.132 


iS) 

\ J 


3.115 


(S) 


3.109 


iS) 


2.060 


^(25) 


3.686 


3.676 


(S) 


3.662 


{S) 


3.635 




3.681 




6.897 


ibilD) 


3.770 


3.794 


(D) 


3.770 


(D) 


3.747 




3.788 




2.060 


V'(35) 


4.040 


4.079 


(S) 


4.049 


(S) 


4.009 


{S) 


4.107 


{S) 


12.636 


ij{2D) 


4.160 


4.156 


{D) 


4.121 


(D) 


4.087 


(D) 


4.168 


(D) 


6.897 


^(45) 


4.415 


4.434 


iS) 


4.386 


(S) 


4.290 


(S) 


4.498 


(S) 


17.227 


^(3D) 




4.482 


{D) 


4.432 


(D) 


4.390 


(D) 


4.581 


(D) 


12.636 






4.781 


(S) 


4.737 


iS) 


4.566 


(S) 


4.812 


(S) 


32.968 






4.889 


(D) 


4.801 


(D) 


4.711 


(D) 


5.181 


(D) 


17.227 






5.135 


iS) 


5.057 


(S) 


4.993 


(S) 


5.429 


(S) 


23.372 






5.451 


(D) 


5.325 


(D) 


5.136 




5.985 


(D) 


32.968 


ij{6D) 




6.030 


(D) 


5.869 


(D) 


5.819 


(D) 


6.943 


(D) 


23.372 



(3) 



Bold numbers stand for the masses included in the fit as an input. The states 1 arc the 
mixture of S and D waves (in parantheses the dominant wave is shown, see indices (nS) and 
{riD)). In Appendix, we present the wave functions for the 1 states and the values Woo, 
Wo2, VF22, which characterize a percentage of the L = and L — 2 components (see [1] for the 
details). The last column gives us the mean square radii for the states in solution 7(cc): i?f 
GeV-2. 

For the other considered states, the fit resulted in the following masses and Rj (all values 
in GcV units): 

for 0"+ states (L = 0): 



State 


Data 


I{cc) 


II{cc) 


U{cc) 


R{cc) 


Rj 


VcilS) 


2.979 


2.979 


2.985 


3.016 


2.970 


1.682 


Vc{2S) 


3.594 


3.606 


3.581 


3.574 


3.616 


6.207 






4.030 


3.994 


3.958 


4.067 


11.813 


Vci^S) 




4.386 


4.338 


4.265 


4.461 


16.604 






4.763 


4.715 


4.555 


4.813 


30.919 






5.136 


5.054 


4.881 


5.475 


22.831 



(4) 



for 0++ states {L = 1): 



State 


Data 


J(cc) 


//(cc) 


U{cc) 


P(cc) 


Rj 


Xco(lP) 


3.415 


3.412 


3.372 


3.473 


3.407 


3.401 


Xco{2P) 




3.867 


3.822 


3.850 


3.876 


8.777 


Xco(3P) 




4.228 


4.180 


4.173 


4.265 


15.115 


Xco(4P) 




4.538 


4.489 


4.493 


4.603 


22.156 


Xco(5P) 




4.943 


4.862 


4.795 


5.038 


18.133 


Xco(6P) 




5.436 


5.318 


5.067 


5.804 


13.806 



(5) 
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for 1++ states = 1): 



State 


Data 


I(cc\ 


Il(rr) 

J. J. \ L-L- 1 


U(rr) 


B(rr) 

-£ LI L-L- 1 


P? 


Xcl(lP) 


3.510 


3.500 


3.491 


3.503 


3.497 


4.234 


Xci(2P) 


3.872 


3.933 


3.904 


3.880 


3.944 


9.861 


Xci(3P) 




4.278 


4.242 


3.989 


4.328 


17.628 


Xci(4P) 




4.563 


4.524 


4.228 


4.809 


24.460 


Xci(5P) 




4.916 


4.841 


4.575 


5.635 


18.407 


Xci(6P) 




5.460 


5.345 


4.819 


6.445 


13.345 



for 1+ states {L = 1): 



State 


Data 


/(cc) 


II{cc) 


U{cc) 


P(cc) 


Rj 


hc{lP) 


3.526 


3.511 


3.513 


3.522 


3.515 


4.447 


hc{2P) 




3.949 


3.918 


4.013 


3.942 


10.199 


hcisP) 




4.323 


4.260 


4.385 


4.300 


14.886 






4.678 


4.599 


4.696 


4.659 


19.976 


K{5P) 




5.158 


4.916 


5.078 


4.997 


24.106 


hc{QP) 




5.968 


5.353 


5.531 


5.472 


15.336 



and for 2++ states (L = 1, 3): 



State 


Data 


/(cc) 




//(cc 




U{cc) 




R{cc) 




Rj 


Xc2(lP) 


3.556 


3.552 


(P) 


3 


570 


in 


3.508 


(P) 


3.556 


(P) 


5.008 


Xc2(2P) 


3.941 


3.986 


(P) 


3 


957 




3.898 


(P) 


3.969 


(P) 


11.085 


Xc2(lP) 




4.050 


(P) 


4 


035 


(F) 


3.946 


(F) 


4.054 


(F) 


5.008 


Xc2(3P) 




4.350 


(P) 


4 


288 


(p) 


4.222 


(P) 


4.320 


(P) 


14.928 


Xc2(2P) 




4.476 


(P) 


4 


333 


(F) 


4.260 


(F) 


4.362 


(F) 


11.085 


Xc2(4P) 




4.786 


(P) 


4 


402 


(P) 


4.546 


(P) 


4.411 


(P) 


41.793 


Xc2(3P) 




5.073 


(P) 


4 


648 


(F) 


4.558 


(F) 


4.705 


(F) 


14.928 


Xc2(5P) 




5.505 


(P) 


4 


687 


(P) 


4.803 


(P) 


4.762 


(P) 


12.018 


Xc2(4P) 




5.833 


(P) 


5 


055 


(F) 


4.937 


(F) 


5.153 


(F) 


41.793 


Xc2(6P) 




6.318 


(P) 


5 


223 


(P) 


5.079 


(P) 


5.334 


(P) 


10.590 


Xc2(5P) 




6.765 


(P) 


5 


596 


(F) 


5.429 


(F) 


5.727 


(F) 


12.018 


Xc2(6F) 




7.710 


(P) 


6 


154 


(F) 


6.065 


(F) 


6.325 


(F) 


10.590 



for 2"+ states (P = 2): 



State Data /(cc) //(cc) U{cc) R{cc) Rj 

r]c2{lD) — 3.829 3.823 3.742 3.818 7.721 

r7c2(2P) — 4.185 4.164 4.087 4.185 14.387 

r/c2(3P) — 4.471 4.449 4.397 4.591 22.729 

7;c2(4P) — 4.788 4.745 4.713 5.151 18.708 

r]c2{^D) — 5.242 5.157 5.084 5.996 14.024 

77c2(6P>) — 5.813 5.693 5.546 6.832 12.227 . 



M(GeVA 



4.4 



4.2 



4.a 



3.8 



3.6 



3.4 



3.2 



3.0 



4P 



4S 



3D 
4S 

2D 



4P 



3P 



2F 
3P 



3P 



4P 



3P 



3D 



2D 




77.(0-^) vcr-) Xco(o^O xc,(i^^) Xo2{2n Vc2i2-1 



Figure 1: The cc levels in solution I{cc). Solid lines stand for observed states, dashed lines for 
the predicted ones. Thin solid lines show the transitions included into fitting procedure, thin 
dashed line demonstrates the transitions, whose widths are predicted. 



6 



In Fig. 1, the levels for solution /(cc) are shown for the mass region M < 4.5 GeV. The 
wave functions for this solution are given in Appendix. 

The obtained variants of solutions allow us to consider the status of X(3872) and y(3941). 

The solution J/(cc) provides for Xci(2P) the mass 3904 MeV that is close to the value 
found in [9, 10, 11, 12] for the pick denoted as X(3872). In agreement with //(cc)-solution, the 
analysis [26] favours the quantum numbers J^'~^ = l"^"*" for this state. 

Still, one cannot exclude that X(3872) is the 2~"'"-state (see [27] and references therein). In 
all solutions, the mass of ?7c2(l-D) is lying in the interval 3815 — 3830 MeV that demonstrates 
the plausibility of this variant as well. 

The signal y(3941) does not contradict the hypothesis about its 2++ nature, because the 
mass of the X2c(2-P) is close to y(3941) in all solutions. 



2.2 Radiative transitions {cc)in ^ 7 + (cc) 



out 



In Fig. 1, we show radiative decays which have been accounted for in the fitting procedure 
(corresonding formulae are presented in [1, 28]). For the levels below DD threshold, experi- 
mental data [25, 29, 30, 31] and the magnitudes of widths obtained in different versions of the 
fit are as follows (in keV): 



Process 



Xco(lP) 

Xcl(lP) 
Xc2(lP) 

r/co(25) ■ 
^(2^) - 
ij{2S) - 
i^{2S) - 
^(25) - 
V^(25) - 



7^co(l^) 
lXci{lP) 



Data 

1.1±0.3 

165±50 

295±90 

390±120 

0.8±0.2 
26±4 
25±4 
20±4 



/(cc) 

2.877 

180.845 

291.536 

226.700 

23.297 

1.661 

22.854 

45.291 

26.096 

0.541 



//(cc) 

4.853 

164.884 

297.362 

239.835 

23.326 

2.529 

32.408 

51.853 

20.929 

0.770 



U{cc) 

1.4 

273.8 

391.8 

312.3 

40.263 

0.37 

12.2 

31.1 

40.2 

1.003 



R{cc) 

4.216 

180.578 

289.569 

225.265 

25.852 

2.309 

22.169 

44.343 

26.620 

0.447 



(10) 



Note that the 20% accuracy is allowed for the transitions ■0(25') — > 7Xcj(l-P) and 30% one for 
Xcj(l-P) 7^(1'5') (we use a bit larger errors than those obtained in [31]). We also predict the 
widths of the decays rjco — > jJ/ip and ip{2S) — > jr]co{2S). 

The calculated values in (10) agree rather reasonably with the data. 

The predictions of widths of the levels above the DD threshold (see Fig. 1) are as follows 
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(widths are in keV): 



Process Data lice) II{cc) U{cc) R{cc) 

Xcoi2P) -fj/ip — 9.181 5.841 0.468 12.612 

Xci(2P)^7J/^ — 48.130 50.831 28.797 59.013 

Xc2(2P) ^ 7J/^ — 56.911 72.080 31.331 67.713 

Xco(2P) ^ 7^(25) — 81.851 73.158 92.450 73.925 

Xci(2P) ^ jip{2S) — 212.937 221.604 290.379 192.502 

Xc2(2P) ^ -fip{2S) — 154.796 176.599 197.162 144.714 



(11) 



2.3 The cc component of the photon wave function and two-photon 
radiative decays 

In the fitting procedure, we approximate the vertex of the transition 7 ^ cc by the following 
formula: 

^ 1 

G,^cc-(5)(«) = E t^n5Gv^(.5)(«) + l + exp(-/?,(. - .0)) ' ^^^^ 
2 

n=l 

where Gv{nS){s) is the vertex for the transition ipijiS) — > cc and Gv(nD){s) is the vertex for the 
transition il){nD) cc, see [1] for the details. The parameters CnSi CnD, l3j, sq for solutions 
/(cc), II{cc), U{cc) and R{cc) have been found as follows (in GeV): 



/(cc 


) 


//(cc) 


U{cc) 


R{cc) 


Cis 


= -10.100 


Cis 


= -9.073 


Cis 


= -3.852 


Cis 


= -14.990 


C2S 


= 2.121 


C2S 


= -0.252 


C2S 


= 0.476 


C2S 


= 3.416 


Css 


= -0.373 


Css 


= 0.038 


Css 


= 0.325 


C3S 


= -0.862 


Cas 


= -1.818 


C45 


= 0.217 


C^s 


= 0.667 


Cas 


= -1.834 


C5S 


= 9.998 


C55 


= -1.158 


C5S 


= -2.571 


C5S 


= 11.704 


Ces 


= 3.369 


Ces 


= -0.529 


Cqs 


= -0.707 


Ces 


= 0.305 


CiD 


= -0.031 


CiD 


= -0.005 


CiD 


= 0.080 


CiD 


= -0.186 


C2D 


= -0.110 


C2D 


= 0.020 


C2D 


= -0.082 


C2D 


= 0.114 




2.85 


6-y = 


: 2.85 


= 


: 2.85 


= 


: 2.85 


So = 


18.79 


So = 


: 18.79 


So = 


: 18.79 


So = 


: 18.79 



(13) 



The corresponding vertices G^^cc{s) are shown in Fig. 2. 

Experimental values of partial widths [25, 32, 33, 34, 35] together with the widths obtained 
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Figure 2: a) Vertices G^^Xcc (a) and Cij^cc (b) for solutions I{cc) (solid lines), U{cc) (dashed 
lines) and R{cc) (dotted lines). 



9 



in the fitting procedure are shown below (in keV): 



Process 





—>■ e^e 


^{2S) - 


^ e+e" 




> e+e~ 


^(35) - 


> e+e^ 


^l^{2D) - 


^ e^e~ 


i^i^S) - 


> e+e~ 



Data 
5.40 ± 0.22 
2.14 ± 0.21 
0.24 ± 0.05 
0.75 ± 0.15 
0.47 ± 0.10 
0.77 ± 0.23 



I{cc) 


II{cc) 


5.392 


5.360 


2.146 


2.137 


0.239 


0.730 


0.754 


0.472 


0.472 


0.236 


0.771 


0.763 



U{cc) R{cc) 

5.403 6.088 

2.142 2.843 

0.240 0.301 

0.749 1.210 

0.469 0.113 

0.770 0.915 



With the determined vertices for G^^cc{s), we can obtain the widths of the two-photon 
decays (see [1] for more detail). The comparison of experimentally measured widths with those 
obtain in our calculations is given below: 



Process 
r/co(15) -f-f 
Xco(l^) ^ 77 
Xc2(l^) ^ 77 



Data -^(cc) II{cc) 

7.0±0.9 6.971 6.808 

2.6±0.5 2.426 2.620 

1.02±0.40±0.17(L3 ) 0.558 0.528 
1.76±0.47±0.40(OPAL) 
1.08±0.30±0.26(CLEO) 
0.33±0.08±0.06(E760) 



U{cc) R{cc) 

7.002 7.088 

2.578 2.212 

0.068 5.706 



(15) 



Let us emphasize that the data do not tell us definitely about the width Xc2(3556) 77. In 
the reaction pp — 77, the value r(x2(3556) — > 77) = 0.32 ± 0.080 ± 0.055 keV was obtained 
in [35], while in direct measurements such as e'^e~ annihilation the width is much larger: 
1.02 ± 0.40 ± 0.17 keV [32] , 1.76 ± 0.47 ± 0.40 keV [33] , 1.08 ± 0.30 ± 0.26 keV [34] . The 
compilation [25] provides us with the value close to that of [35]. The value found in our fit 
agrees with data reported by [32, 33, 34] and contradicts the magnitude from [35]. 

Our predictions of widths cc — > 77 for the levels below 4 GeV are as follows: 

Process Data I{cc) II{cc) U{cc) R{cc) 

r]coi2S) ^-f-f — 12.234 11.947 12.289 12.439 

Xco(2P)^ 77 — 2.142 2.313 2.276 1.953 (16) 

Xc2(2P)^ 77 — 0.500 0.473 0.061 5.111 



In tables 1,2,3, we compare our results with those obtained by other authors. 

In [36], the calculated widths depend on a chosen gauge for the gluon exchange interaction 
— we demonstrate the results obtained for both Feynman (F) and Coulomb (C) gauges. 

In [37], the cc system was studied in terms of scalar (S) and vector (V) confinement forces — 
both variants are presented in tables 1,2,3. The results obtained in the nonrelativistic approach 
to the cc system [38] are also shown in tables 1,2,3. 
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In both rclativistic [36, 37] and nonrelativistic [38] approaches, there is rather large dis- 
crepancy between the data and calculated values of ip{nS) e+e^ (in [36] the width of the 
transition J/ip ^ e'^e" was fixed with the use of a subtraction parameter). In our opinion, the 
reason is that in all above-mentioned papers, soft interaction of quarks was not accounted for 
— we mean the processes shown in Fig. 3b, c of Ref. [1]. In fact, the necessity of taking into 
consideration the low-energy quark interaction was imderstood decades ago but untill now this 
procedure has not become commonly accepted even for light quarks; see, for example, [39, 40]. 

Table 1: Comparison of data for the transitions (cc)j„ ^7-1- {cc)out with our results and 
calculations of other groups (the width is given in keV). 



DccaA- 


Data 


lice) 


LS(F)[36] 


LS(C)[36] 


KM(R)[37] 


KM(V)r37] 


>JK[.38] 




- '/cu(l.S')- 


i.iio.:-! 


:-!. i 


1.7 1.:-! 


1.7 1.1 


:13.5 


2.(>(i 


1.21 


V(25) ^ 


Xco(lP)7 


26±4 


23 


31-47 


26-31 


31 


32 


19.4 


^P{2S) ^ 


Xcl(lP)7 


25±4 


45 


58-49 


63-50 


36 


48 


34.8 


^{2S) ^ 


Xc2(lP)7 


20±4 


26 


48-47 


51-49 


60 


35 


29.3 


V(25) ^ 


Vco{lS)-i 


0.8±0.2 


2.0 


11-10 


10-13 


6 


1.3 


4.47 


Xco(lP) 


^ J/V(15')7 


165±50 


181 


130-96 


143-110 


140 


119 


147 


Xcl(lP) 


^ J/V(1S)7 


295±90 


292 


390-399 


426-434 


250 


230 


287 


Xc2(lP) 


j/V{is)t 


390±120 


227 


218-195 


240-218 


270 


347 


393 



Table 2: Decay widths e+e , our results and those of other groups v.s. data (in keV). 



Decay 


Data 


I(cc) 


LS(F)[36] 


LS(C)[36] 


RM(S)[37] 


RM(V)[37] 


NR[38] 


J/^{1S) 




e+e 


5.40 ± 0.22 


5.39 


5.26 


5.26 


8.05 


9.21 


12.2 


V(25) ^ 


e^ 




2.14 ± 0.21 


2.15 


2.8-2.5 


2.9-2.7 


4.30 


5.87 


4.63 


V.(1D) ^ 


e 


<-e- 


0.24 ± 0.05 


0.24 


2.0-1.6 


2.1-1.8 


3.05 


4.81 


3.20 


V(35) ^ 


e^ 


"e^ 


0.75 ± 0.15 


0.75 


1.4-1.0 


1.6-1.3 


2.16 


3.95 


2.41 


V>(2D) ^ 


e 


^e- 


0.47 ± 0.10 


0.47 












^-(45) ^ 


e^ 


"e~ 


0.77 ± 0.23 


0.77 













Table 3: Decay widhs (cc)i„ — > 7 -|- {cc)out in the region below 4 GeV (in keV) 



Dctay 


Data 


l[cc) 




lllj 


M 


113| 






116] 


r)c{lS) 77 


7.0±0.9 


6.9 


6.2-6.3 (F,C) 


5.5 


3.5 


10.9 


7.8 


5.5 


4.8 


r)c{lS) 77 




12.2 




1.8 


1.38 




3.5 


2.1 


3.7 


Xco(l-P) 77 


2.6±0.5 


2.4 


1.5-1.8 (F,C) 


2.9 


1.39 


6.4 


2.5 


5.32 




Xco(l-P) 77 




2.1 




1.9 


1.11 










Xc2(lP) -> 77 


1.02±0.40±0.17[32] 


0.61 


0.3-0.4 (F,C) 


0.50 


0.44 


0.57 


0.28 


0.44 






1.76±0.47±0.40[33] 




















1.08±0.30±0.26[34] 




















0.33±0.08±0.06[35] 


















Xc2(lf ) -> 77 




0.5 




0.52 


0.48 
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3 Conclusion 



We have performed a successful description of both cc levels and their radial excitation tran- 
sitions for several interaction variants as follows: for scalar confinement potential, scalar and 
vector confinement potential in case of instantaneous and retardation forces. Such a diversity 
in the choice of forces reflects the lack of available experimental data to guarantee unambigu- 
ous determination of the interaction. First of all, one misses the data for the radiative cc 
decays. This is why we pay special attention to the predictions given by different versions of 
our calculations. 

Rather good description of the cc levels being obtained, the wave functions related to dif- 
ferent solutions may noticeably differ from each other, and sometimes they demonstrate rather 
unusual behaviour. Because of that, we pay a considerable attention to the presentation of wave 

functions of the cc systems. Wc think that the future progress in understanding the c-quark 
interactions in the soft region is related precisely to a better knowledge of wave functions, that 
should be based on further study of radiative transitions. 
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4 Appendix. Wave function of the cc sector 

** 

Tables 4-10 give us the Ci{S, L, J; n) coefficients, which determine the wave functions 
for solution /(cc) according to the following formula: 

i=l 

where = (recall that s = Arn? + 4k^). The fitting parameter j3 is of the order of 0.5 — 1.5 
GeV^^ and it may be different in different flavor sectors; we put /? = 1.2 GeV^^. In tables, we 
also show Wll''s, which enter the normahzation condition Wll + Wll> + Wl'l> — 1, being the 
convolution of the wave functions 

see [1] for more detail. In Figs. 3-9, we demonstrate these wave functions. 
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Table 4: Constants Ci{S, L, J;n) from Eq. (17) (in GeV units) for the wave functions of T- 
mesons in solution /(cc) 



^(15) 



IV'oo Wo2 

0.1)9987 -0.0()0()(i 



\V>> 
0.00079 



^(25) 



H Of) W()2 
0.99921 -0.00171 



0.002."')0 



Won 
0.99807 



-0.00288 



11 22 

0.00181 



7.9009 
-6.0655 
16.7028 

-105.4498 
224.9337 

-231.2288 
127.1267 
-36.0487 
4.1617 



-5.9190 
34.4862 
-83.5239 
102.1413 
-63.9523 
16.2702 
1.9716 
-1.8633 
0.2684 



11.6473 
47.6240 
-453.2935 
1094.1453 
-1334.8514 
943.6424 
-393.8494 
90.4811 
-8.8655 



1.5241 
-14.8598 
41.5693 

-41.5148 
2.7143 
24.1754 
-18.1178 
5.3507 
-0.5785 



17.3698 
41.6854 
-860.2544 

2693.3655 
-3793.6194 
2870.7394 
-1206.4102 
264.5053 
-23.4715 



-30.6053 
199.5601 
-499.8739 

625.7537 
-415.9994 
134.3492 
-9.1646 
-5.3519 
0.9886 



^(45) 



*(55) 



*{6S) 



Woo 
0.99166 



Wo2 

-0.00737 



W22 
0.01570 



Woo 
0.99889 



W02 
-0.00450 



W22 
0.00561 



Woo 
0.95583 



W02 
-0.02419 



W22 
0.06835 



36.5862 
-212.7589 
-6.8022 
1855.8350 
-4206.0402 
4182.8669 
-2151.4758 
559.1998 
-58.1665 



-1.1743 
51.7345 
-243.2839 
448.2876 
-401.5179 
178.7153 
-31.8182 
-1.4693 
0.8430 



-95, 
1136. 
-4872. 
10215. 
-11814. 
7907. 
-3052. 
631. 
-54. 



0826 
7086 
0833 
8015 
1052 
6179 
3529 
1276 
2492 



-9.0172 
81.4929 
-284.5462 
523.5176 
-570.3848 
381.1356 
-152.6765 
33.4538 
-3.0632 



78, 
-1113, 
5651, 
-13997, 
18993, 
-14736, 
6486, 
-1501, 
141, 



7073 
5476 
5790 
7499 
0373 
1180 
3763 
0222 
2673 



62.0712 
-481.43.54 
1419.5113 
-2062.1487 
1573.0804 
-594.8321 
71.5608 
15.2345 
-3.6824 



^(ID) 



^(2i3) 



^(3D) 



Woo 
0.00435 



W02 
-0.00169 



W22 
0.99734 



Woo 
0.00736 



W02 
-0.00316 



W22 
0.99581 



Woo 
0.00697 



W02 
-0.00433 



W22 
0.99736 



^(1,0,1) 



1,2,1) 



1,2,1) 



84.4242 
-205.1653 
181.7717 
-77.1585 
49.8493 
-57.6810 
35.7092 
-10.1511 
1.0906 



-0.1030 
3.1892 
-19.8052 
46.4276 
-58.9700 
44.8041 
-20.4094 
5.1363 
-0.5496 



313.1175 
-1319.9899 
2113.8219 
-1592.4362 
460.0515 
111.8270 
-120.5831 
32.7133 
-3.1245 



-0.2904 
10.2284 
-70.5461 
192.4214 
-260.8978 
194.4329 
-80.7753 
17.4329 
-1.5051 



467.5023 
-2389.2054 
4394.0661 
-3418.8591 
524.5556 
907.2750 
-629.5731 
163.7261 
-15.6276 



-1.8858 
23.4791 
-121.9462 
332.6430 
-507.1882 
442.0815 
-218.2705 
56.6805 
-6.0100 



^(4D) 



*(5D) 



^(6D) 



Woo 
0.02465 



W02 
-0.01134 



W22 
0.98669 



Woo 
0.01698 



W02 
-0.01676 



W22 
0.99978 



Woo 
0.02212 



W02 
-0.01885 



W22 
0.99673 



^(1,0,1) 



413.1586 
-2384.9885 
4895.6915 
-4120.8799 
490.3119 
1523.0854 
-1076.5260 
293.9187 
-29.5824 



8.8245 
-114.1754 
524.1161 
-1159.3614 
1390.2561 
-947.6269 
366.2696 
-74.6692 
6.2449 



-289.4977 
1734.9402 
-3619.8761 
2815.1704 
415.7821 
-2036.6024 
1288.9620 
-343.3866 
34.1359 



1.5578 
-22.5966 
127.1963 
-371.5994 
619.5920 
-602.2710 
332.9247 
-96.0393 
11.1873 



244.8123 
-1545.0640 
3437.6769 
-2942.2393 
-231.3078 
2146.9545 
-1496.2120 
431.2140 
-46.2070 



3.7011 
-63.3077 
378.8322 
-1087.4587 
1690.6255 
-1488.8664 
736.8267 
-189.8749 
19.7831 
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Table 5: Constants Ci{S,L, J;n) from Eq. (17) (in GeV) for the wave functions of Xc2-inesons 
in solution /(cc) 





Xc2(lP) 


Xc2(2P) 


Xc:2 


,3P) 






W'li H 




Ws.i 


Wri 11 


'1.3 113.3 


11 '11 11 


1.3 


ll:i3 




0.1)91)52 -0.00070 


0.00118 


0.99705 -0.00211 0.00 IKi 


0.99209 -0.00518 


0.01:«)9 




y,(i,l,:^) 


^(1,3,2) 




^(1,3,2) 


^(1,1,2) 


^(1,3,2) 




-17.5035 




9.0856 


4.5957 


39.0522 


-93.1247 




44.4200 


2 


-86.9407 




-40.9452 


-623.9524 


-250.9432 


163.8698 




-311.5971 


3 


506.5360 




85.1047 


3316.2121 


629.3054 


1317.7565 




796.6134 


4 


-988.9820 




-100.6670 


-7222.9252 


-821.2865 


-5066.0027 




-990.8076 


5 


1026.3189 




71.0253 


8391.4505 


608.6930 


7472.6180 




644.8615 


g 


-625.3392 




-29.4658 


-5639.3954 


-257.5492 


-5744.7189 




-202.3282 


7 


224.1992 




6.6484 


2197.9985 


57.7459 


2434.8192 




13.3976 


8 


-43.7175 




-0.6473 


-460.8963 


-5.2652 


-539.0269 




7.5182 


9 


3.5664 




0.0062 


40.1760 


-0.0143 


48.6569 




-1.3447 




Xc2{4P) 


Xc2(5P) 


Xc2(6P) 




Wii Wis 


Wss 


Wu Wis Wss 


Wll Wis 


W33 




0.88640 -0.00477 


0.11837 


0.98493 -0.01273 0.02780 


0.93818 -0.03946 


0.10128 


2 


^■(1,1-2) 


^(1.3.2) 


^^(1.1,2) 


^(1.3,2) 


.^,(1.1.2) 


^(1,3,2) 


1 


638.6609 




391.6378 


-134.7570 


-221.6602 


117.7874 




104.0293 


2 


-5864.8034 


-2008.6939 


1751.9157 


1250.1138 


-1488.3254 




-706.9778 


3 


21108.8607 


4224.3629 


-8121.0456 


-2766.6069 


6956.8543 


1877.3877 


4 


-39600.9145 


-4709.0719 


18375.8266 


3074.1617 


-16174.3244 


-2509.6472 


5 


42944.0995 


2957.4768 


-22896.3579 


-1771.0161 


20877.3955 


1793.4813 


Q 


-27861.8978 




-999.8571 


16466.2322 


434.8957 


-15563.0676 




-647.3392 


7 


10649.1787 




137.5405 


-6796.9738 


31.0544 


6631.3660 




78.8171 


8 


-2206.5653 




8.3919 


1492.7938 


-36.5787 


-1494.2661 




13.2404 


9 


190.7603 




-3.1287 


-134.9738 


5.1157 


137.6804 




-3.2854 




Xc2(lF) 


Xc2(2P) 


Xc2(3P) 




Wu Wis 


Wss 


Wll Wis Wss 


Wu Wis 


Wss 




0.00301 -0.00135 


0.99834 


0.10992 -0.00454 0.89462 


0.03903 -0.01166 


0.97263 


2 


^(1.3.2) 


^(1,1,2) 


,^(1.3.2) 


^(1.1,2) 


^(1,3,2) 


^(1,1.2) 


I 


-242.7471 




10.5652 


-565.8475 


225.5515 


733.1238 




19.8324 


2 


744.2459 




-95.8214 


2117.9502 


-2067.2776 


-3370.3257 




-120.2819 


3 


-845.6564 




355.4159 


-2625.9078 


7468.3409 


5512.4695 




179.6152 


4 


320.2998 




-686.5646 


638.7903 


-14099.5057 


-3448.4647 




157.0118 


5 


166.8886 




763.7615 


1409.3867 


15396.8275 


-410.2084 




-688.2629 


g 


-224.3091 




-506.8350 


-1486.9023 


-10057.8914 


1736.5500 




755.6477 


7 


96.3257 




197.7245 


638.6645 


3868.9470 


-956.7876 




-394.8955 


8 


-19.3136 




-41.7533 


-131.7223 


-806.4039 


226.5196 




101.4201 


9 


1.5203 




3.6750 


10.7017 


70.0913 


-20.3853 




-10.2876 




Xc2(4F) 


Xc2(5F) 


Xc2(6P) 




Wll Wl3 




Wll Wr.i VK33 


Wll Wl3 


1^33 




0.01125 -0.00896 


0.99771 


0.01670 -0.01856 1.00186 


0.01548 -0.02033 


1.00485 


i 




^,(1.1.2) 


^(1,3,2) 


^(1.1,2) 


^(1,3,2) 


V,(l.l>2) 


1 


-580.3877 




-12.4060 


431.6178 


1.5693 


345.9078 




5.2688 


2 


2948.0418 




137.4061 


-2289.0231 


-16.1472 


-1924.9704 




-74.4252 


3 


-5368.6253 




-566.3875 


4363.4753 


50.9676 


3911.9978 




376.9102 


4 


3893.8813 


1161.4946 


-3306.3734 


-47.7121 


-3311.2035 




-935.6268 


5 


3.8221 


-1320.6326 


-27.4883 


-46.6674 


358.9932 


1279.4826 


6 


-1709.0721 




866.8103 


1606.1715 


127.5496 


1343.8856 


-1004.6984 


7 


1047.0953 




-325.6228 


-1012.8337 


-98.5000 


-964.7262 




448.1251 


8 


-263.5517 




64.7935 


261.2351 


33.0133 


268.3202 




-104.9359 


9 


24.8409 




-5.2798 


-25.0732 


-4.1196 


-27.3899 




9.9921 
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Table 6: Constants Ci{S, L, J; n) from eq. (17) (in GeV) for the wave functions of ?7co, Xco, Xci 
and hci mesons in solution /(cc) 



i 




?7co(2P) 


'?co(3P) 


r?co(4P) 


??co(5P) 


??co(6P) 


1 


7.6494 


-10.5636 


18.1089 


33.6294 


-89.2312 


-75.5014 


2 


-23.0211 


-53.9217 


15.0118 


-170.0167 


1032.3991 


1029.1523 


3 


124.2754 


461.0285 


-623.7988 


-168.7223 


-4276.2522 


-5025.8538 


4 


-395.6413 


-1133.8010 


1908.5781 


2060.7362 


8662.4563 


11959.9547 


5 


656.5275 


1456.7298 


-2536.9289 


-4211.5707 


-9688.9143 


-15569.2648 


6 


-605.5583 


-1101.1359 


1780.5781 


3997.5023 


6289.9521 


11572.6965 


7 


314.7296 


493.4357 


-680.8186 


-1997.5324 


-2365.3552 


-4872.6646 


8 


-86.2318 


-121.2544 


132.0261 


509.0043 


479.4585 


1076.4010 


9 


9.7057 


12.6004 


-9.8555 


-52.2319 


-40.7532 


-96.3877 


i 


Xco(l-P) 


Xco(2P) 


Xco(3P) 


Xco(4P) 


Xco(5P) 


Xco(6P) 


1 


-17.8495 


-50.2837 


-151.3174 


-313.0348 


318.7744 


-249.2578 


2 


3.8696 


5.7940 


731.5190 


2334.3640 


-2879.7647 


2448.4531 


3 


29.1978 


634.0603 


-1074.8942 


-6487.9099 


9875.5585 


-9126.1510 


4 


70.7786 


-1637.8698 


270.8311 


8844.2446 


-17101.9307 


17108.1285 


5 


-273.1091 


1949.8183 


727.7711 


-6388.0162 


16610.2166 


-17802.3919 


6 


321.4758 


-1315.0808 


-767.5997 


2335.4592 


-9396.2245 


10577.7296 


7 


-184.1650 


521.0477 


313.2245 


-289.9459 


3061.3295 


-3502.4000 


8 


52.5331 


-113.9435 


-55.4903 


-50.4060 


-530.9090 


584.2914 


9 


-5.9970 


10.6960 


3.0610 


13.0531 


37.9718 


-36.1088 


i 


Xcl(lP) 


Xci(2P) 


Xci(3P) 


Xci(4P) 


Xci(5P) 


Xci(6P) 


1 


23.0070 


56.7513 


217.3310 


-388.4208 


344.6565 


238.1913 


2 


-3.8670 


-14.0364 


-1332.2767 


3249.2463 


-3365.3754 


-2490.3066 


3 


-115.4587 


-778.7437 


3058.3410 


-10453.2199 


12574.5044 


9950.4054 


4 


191.1470 


2115.0575 


-3497.5882 


17295.2284 


-23976.3393 


-20229.2674 


5 


-112.1101 


-2567.5452 


2172.1531 


-16412.7705 


25956.4509 


23216.4628 


6 


4.4321 


1719.1521 


-716.3965 


9284.6416 


-16581.9431 


-15583.6367 


7 


23.5865 


-657.5312 


101.6710 


-3089.5089 


6182.0939 


6036.8548 


8 


-9.8945 


134.9058 


1.4698 


556.8344 


-1242.3831 


-1245.5237 


9 


1.2824 


-11.5449 


-1.2720 


-41.8049 


103.9046 


105.5852 


i 






hAm 


/'.:(4P) 


/',:(5P) 


//,;(()P) 


1 


-19.4824 


44.2886 


130.8468 


-296.9419 


431.5082 


-293.9990 


2 


-42.1580 


119.7469 


-553.3042 


2284.4817 


-4002.0767 


2983.8574 


3 


305.3320 


-1319.3114 


298.1904 


-6661.4528 


14363.1096 


-11699.4557 


4 


-573.6405 


3202.4192 


1602.7564 


9794.3567 


-26533.5962 


23525.9355 


5 


548.0217 


-3792.3619 


-3269.2777 


-8048.5695 


28014.7411 


-26866.9716 


6 


-299.4493 


2527.4497 


2752.3010 


3805.6849 


-17540.3649 


18029.3913 


7 


93.8170 


-967.1883 


-1199.8706 


-1003.0218 


6432.1760 


-7008.3077 


8 


-15.4547 


198.4952 


266.1679 


130.1848 


-1274.7649 


1455.2756 


9 


1.0081 


-1().9394 


-23.7422 


-3.G5G9 


105.3324 


-124.497() 
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Figure 3: Wave functions for ip{nS) in solution /(cc). Solid lines stand for V'^^'"'^-* and dashed 
ones for 'ip^''^'^\ All magnitudes are in GeV units. 
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Figure 4: Wave functions for ip{nD) in solution /(cc). Solid lines stand for V'^^'^'^-* and dashed 
ones for All magnitudes are in GeV units. 
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Figure 5: Wave functions for r]c in solution 7(cc). 
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Figure 6: Wave functions for Xco in solution 7(cc). 
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Figure 7: Wave functions for Xci in solution 7(cc). 
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Figure 8: Wave functions for Xc2 in solution /(cc). Solid lines stand for ^^^'^'"^^ and dashed ones 
for ^(i'3'2). 
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Figure 9: Wave functions for Xc2 in Solution /(cc). Solid lines stand for ipl;^''^'"^^ and dashed ones 
for ^^I'l'^l 
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